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IMPLICATIONS  OF  VORTEX  THEORY 
FOR  FIREBALL  MOTION 


1.0  INTRODUCTION 

This  is  an  interim  report  which  collects  the  author's  present 
ideas  on  the  state  of  vortex  theory  and  such  implications  as  can  be  drawn 
regarding  low  altitude  fireball  rise  and  development.  The  contents  of  the 
following  discussion  are  about  half  original  and  half  arguments  from  a 
variety  of  sources  covering  100  years  of  research,  recast  as  to  language 
and  emphasis.  It  is  likely  that  the  reader  will  find  the  discussions  some- 
what confusing  and  unsatisfactory.  That  is  because  neither  the  author  nor 
anyone  else  understands  vortex  motion. 

The  collection  of  this  material  is  believed  to  be  worth  publica- 
tion despite  the  unsatisfactory  status  of  the  theory  because  there  are  a 
few  fairly  firm  facts,  not  all  of  which  are  generally  known,  and  a larger 
number  of  fairly  strong  inferences  which  can  be  drawn  and  which  should  be 
helpful  to  theorists  and  modelers  of  fireball  and  atmospheric  "thermal" 
motion.  It  is  the  author's  hope  that  this  report  can  provide  a helpful 
starting  point  for  future  research  on  low  altitude  fireballs. 


The  fo! lowing  sections  proceed  systematically  from  a general  state- 
ment of  the  low  altitude  fireball  rise  problem  through  a definition  of  the 
concepts  to  be  used,  development  of  some  standard  vortex  theory  results, 
and  a series  of  extensions  of  the  theory  toward  the  case  of  interest.  It 
has  been  necessary  to  make  some  indefensible  assumptions,  relying  on 
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Preceding  page  blank 


intuition  and  reasoning  by  analogy.  It  is  hoped  that  enough  of  the  results 
are  reasonable  to  allow  identification  and  correction  of  the  unreasonable 
ones  in  future  work. 

2.0  GENERAL  CHARACTERISTICS  OF  BUBBLE  RISE  MODELS 

Whenever  the  pressure  equilibrium  radius  of  a fireball  exists 
and  is  small  compared  to  the  atmospheric  scale  height,  it  is  traditional 
to  treat  the  rise  and  stabilization  of  the  fireball  by  a "bubble  rise" 
model . 


A number  of  such  models  have  been  developed  and  can  be  found 
in  the  literature. 


Assumptions  common  to  all  of  the  models  are: 

1.  Bubble  pressure  is  equal  to  local  ambient  pressure. 

2.  The  bubble  entrains  ambient  air  mass  at  a rate 
proportional  to  the  rise  velocity,  some  area 
associated  with  the  bubble,  the  ambient  density, 
and  an  "entrainment  coefficient",  a. 

Usually  the  bubble  is  assumed  to  be  spherical  in  shape. 

Such  models  always  meet  immediately  with  a limited  success,  in 
the  sense  that  a reasonable  value  of  a can  be  found  such  that  stabiliza- 
tion altitude  can  be  predicted  within  a few  kilometers  for  a wide  range  in 
yield  and  some  range  in  burst  altitude. 

Weaknesses  in  such  models  are:  (1)  it  is  not  clear  that  a should 

be  constant  nor  what  parameters  determine  it,  (2)  if  dynamic  informatior. 
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(e.g.,  velocity)  is  required,  another  set  of  rather  arbitrary  assumptions 
is  added,  (3)  factors  believed  to  be  important  have  never  been  treated,  as 
real  atmospheric  temperature,  humidity,  wind  shear  profiles,  turbulence, 
and  the  detailed  profile  of  temperature,  density,  etc. 

The  basic  reason  for  "he  above  situation  is  that  one  is  attempt- 
ing to  extract  desired  results  from  a complex  multidimensional  turbulent 
hydrodynamic  problem  without  recourse  to  a detailed  solution,  and  without 
enough  conservation  laws. 

In  particular,  one  normally  has:  (1)  either  momentum  or  energy 

conservation,  (2)  mass  conservation,  and  (3)  the  condition  of  pressure 
equilibrium  to  determine:  (1)  density,  (2)  temperature,  (3)  volume,  and 

(4)  entrainment  parameter.  Hence  one  is  forced  to  make  an  additional 
assumption.  In  cases  where  both  momentum  and  energy  conservation  are 
used,  a fifth  unknown  appears,  the  velocity,  so  one  is  still  one  equation 
short  and  must  invent  two  arbitrar\  relations,  one  for  a and  one  for  drag 
coefficient,  etc.  rather  than  just  one  (for  a). 

Some  workers  have  attempted  to  resolve  this  dilemma  by  appealing 
to  conservation  of  circulation,  in  particular  the  Mt.  Auburn  group  has  been 
a proponent  of  this  method  of  resolution.  We  will  follow  the'r  lead  but 
with  great  trepidation. 

The  reason  for  our  caution  is  that  almost  all  useful  circulation 
theorems  apply  only  to  incompressible  flow  in  uniform  media,  certainly  not 
a situation  of  interest  to  us.  Circulation  is  defined  as 


r = J (v*v)  • ds  = $ v-d  t 


(1) 


and  accordingly  is  a kinematical  quantity,  not  a dynamical  one  (the  density, 
p,  does  not  appear  in  Equation  1).  If  p is  constant  and  everywhere 
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uniform*  the  distinction  is  academic  and  T is  a useful  quantity,  but  its 
utility  is  suspect  whenever  p is  variable. 

A rationale  exists  which  allows  the  "impulse"  of  a vortex  to  be 
used  as  a physically  meaningful  quantity.1  For  a spheric; 1 vortex  the 
equation  is 


where  P is  impulse,  V is  volume,  U is  mean  vortex  rise  velocity  and 
p is  the  density,  assumed  constant  everywhere.  The  form  Equation  2 
associates  momentum  pVU  with  the  mass  of  fluid  which  accompanies  the 
vortex  in  its  motion  and  takes  account  of  momentum  which  must  be  radiated 
to  infinity,  j pVU,  in  order  to  establish  the  external  flow  field. 

For  the  case  of  a real  atmosphere  the  momentum  and  energy  equa- 
tions become  more  mysterious,  in  that  both  integrals  may  receive  major 
contributions  from  the  air  far  below  the  region  of  interest  due  to  the 
exponential  stratification  of  air  density. 

The  theoretical  status  of  circulation  and  related  quantities 
associated  with  real  vortex  motion  is  thus  in  a very  unsatisfactory  state. 
However,  considerable  research  has  been  done  recently  in  applications  of 
such  concepts  to  atmospheric  thermals  by  dynamical  meteorologists  ’ with 
considerable  success.  It  therefore  appears  worthwhile  to  attempt  to  use 
the  concept  of  circulation  to  furnish  an  extra  relation  which  allows  the 
entrainment  parameter,  a,  to  be  determined  as  a function  of  ether 
bubble  parameters,  including  circulation.  If  successful,  the  determination 


* This  is  assumed  in  almost  all  vortex  theory. 
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of  a then  follows  from  a determination  of  all  of  the  ways  in  which  cir- 
culation can  be  alt  red.  One  hopes  that  these  various  ways  can  be  tabulated 
and  successfully  attacked  one  by  one. 

3.0  DEFINITIONS  OF  CONCEPTS 

3.1  VORTEX 

For  our  purposes  a vortex  is  defined  as  a mass  of  air  which 

contains  vorticity  (Vxv  ^ 0)  and  which  moves  as  a unit  through  the 

in ■ rounding  medium.  We  will  assume  it  to  be  defined  by  a simple  surface, 
as  a spheroid,  whose  position  is  determined  by  the  condition  that  material 
does  not  cross  the  surface.  A sketch  of  a vortex  is  shown  in  Figure  1. 

Note  the  vortex  is  not  identical  to  the  visible  torus,  when 

such  exists.  The  visible  torus  is  included  in  the  vortex.  We  have  shown 
the  surface  as  closed.  For  vortices  which  entrain,  the  surface  is  believed 
to  be  open  at  the  bottom  where  material  enters.  It  is  likely  that  in  many 
cases  the  surface  is  unstable,  resulting  in  entrainment  (and/or  loss)  of 
vortex  material  across  the  boundary.  The  usual  way  of  treating  this  pro- 
cess is  to  continuously  redefine  the  surface  position.3  Thus  the  vortex 
boundary  is  frequently  troublesome  to  define.  This  problem  does  not  affect 
the  concept  of  the  vortex  being  defined  as  all  of  the  material  moving  as  a 
unit  with  the  vorticity,  nor  does  it  reduce  the  value  of  the  concept  of  the 
surface. 

3.2  HILL  VORTEX 

The  only  vortex  form  which  is  stable  and  has  been  described  with 
simple  comnlptc  equations  is  Hill's  spherical  vortex.4  This  vortex  moves 
through  a uniform,  incompressible  medium  with  constant  velocity  and  no 
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entrainment  and  there  are  n£  external  forces  acting.  It  is  far  from  the 
situation  cf  interest  for  fireball  motion*  but  the  only  description  avail- 
able which  is  easy  to  work  with. 


Inside  the  Vortex 


The  vorticity,  w , is  given  by  (see  Figure  1), 

03  = -y-  Ur  sin6/a2  . (3) 

In  the  frame  moving  with  the  vortex  the  velocity  components  are, 
v^  = j U[l-r2/a2 (l+sin20) ] , (4) 

v = j U sin0  cos0  r2/a2  , (5) 

vr  = | U cos0  (l-r2/a2)  , (6) 

and 

vQ  = | U sin0  (2rz/a2-l)  . (7) 

Outside  the  Vortex 

oj  = 0 . (8) 

In  the  frame  moving  with  the  vortex 

v7  = U[l- (1-|  sin20)  a3/r3]  , (9) 

* A count  of  the  underlined  words  in  the  preceeding  sentence  reveals  no 
less  than  five  assumptions  which  are  invalid  for  the  case  of  interest. 
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(10) 


and 


v = 4 U sinQ  cos6  , 
y * r 

= - U cost}  (l-a3/r3)  , 

vQ  = U sin6  (1+j  a3/r3) 


Gross  Quantities 


Mean  square  radius  of  vorticity: 


<r2> 

U) 


Impulse: 

P = 2tt  p a3U 


Energy: 

r-  10  „ 3,,2 

E = — 7r  p a U 
Circulation: 

T = 5 U a 

4.0  DIMENSIONAL  ANALYSIS  FOR  RISING  VORTEX 


We  begin  with  some  fairly  general  results. 

4.1  Entrainment  Parameter  and  Radius 


(ID 

(12) 


(13) 

(14) 

(15) 
(lb) 


Under  any  of  the  three  assumptions:  (a)  incompressible  flow, 

(b)  pressure  equilibrium  without  smali.  scale  mixing,  or  (c)  pressure 
equilibrium  with  vortex  gas  constant,  y , the  same  as  the  exterior  gas 
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constant,  the  increase  of  volume  due  to  entrainment  of  materials  is 
given  by 


dV  = £ dM  (17) 

where  V is  volume,  p&  is  outer  material  density  and  M is  mass  of 
vortex. 


For  a spherical  vortex  we  can  always  write 
dM  5 

^ = 47Ta2p  a (18) 

dz  a v J 

where  z is  altitude,  a is  radius,  and  a is  the  entrainment  parameter, 
defined  by  Equation  18. 

But  by  differentiation  of  the  LHS  of  Equation  17,  it  follows 
that  for  spherical  vortices  and  any  of  assumptions  a-c 


da 

dz 


a 


(19) 


4.2  Constancy  of  Buoyant  Force  in  Uniform  Medium 

The  buoyant  force  operating  on  the  vortex  is  given  by 

Fb  = g (pa-p)  j tt  a3  , (20) 

where  g is  the  acceleration  of  gravity. 
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If  p is  constant  then 
a 


dz 


& . 3a2 
dz 


(pa-p)  at) 


(21) 


but 


dp 

dz 


= 3p  a/a 
a 


p dV 
V dz 


(22) 


where  the  first  term  (derivative  of  M)  follows  from  Equation  18. 
Expanding  the  last  term  in  Equation  22  yields 


dp  _ 
dz 


3p  a/a  - 3pa/a 

3. 


(23) 


Substitution  of  Equation  23  into  Equation  21  shows  the  buoyancy 
force  to  be  constant 


_JU  Ittr  [-  3aa1 2(pa-p)  + 3aa2(Pa-p)]  = 0 (24) 

Almost  every  useful  result  connected  with  vortex  motion  is 
dependent  upon  buoyant  force  being  constant,  most  analyses  in  fact  assume 
zero  buoyancy. 

Here  we  summarize  the  assumptions  which  we  have  made  so  far, 
and  which  underly  most  literature  on  the  subject  of  vortex  motion: 

1.  Time  similarity.  This  assumption  allows  us  to  characterize 
the  vortex  by  a single  average  value  of  such  physical 

parameters  as  radius,  mass,  volume,  force,  etc. 
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2.  One  of  (a)  incompressibility,  (b)  pressure  equilibrium 
without  small  scale  mixing,  (c)  pressure  equilibrium  with 
interior  and  exterior  gas  constants  equal. 

3.  Spherical  shape.  This  is  not  an  important  assumption, 
the  result  (liquation  19)  is  only  slightly  modified  if 
one  allows  a more  general  shape  and  the  result  of  Equa- 
tion 24  is  unchanged.  It  is  an  assumption  of  convenience. 

4.3  Dimensional  Analysis 

We  now  add  one  more  assumption  which  is  almost  universal  in  the 
literature,  and  is  valid  for  most  practical  cases.  We  make  the  approxi- 
mation ; 

4.  That  in  the  equations  of  motion  we  treat  the  vortex  den- 
sity, p,  as  if  it  were  equal  to  the  outer  density,  p , 
in  all  terms  where  it  appears  by  itself,  rather  than  as 

the  difference  p -c  . 

a 

I'or  example,  in  the  impulse  equation  for  a Hill  spherical  vortex, 

j1-  ( 2tt p a 3 U ) = g(Pa-p)  | ira3  , (25) 

we  are  allowed  to  eplace  the  p on  the  left  hand  side  with  p , which 

u 

restricts  us  to  consider  cases  where  the  average  density  inside  the  vortex 
is  only  slightly  different  from  p,  . 

cl 

Having  done  this  we  can  write 

(a 3 U ) = |g  (Pa-p)a3/paEF*,  (26) 


and  by  Equation  24,  the  driving  force,  F*,  is  a constant  with  the 
dimens ions  [L4T  2 ] . 
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In  searching  for  a time  similar  solution  for  U and  a,  we 
can  imagine  only  two  independent  variables  which  might  be  of  physical 
significance,  z,  and  F*  (in  reality,  this  is  a matter  of  taste,  we 
don't  propose  to  defend  this  statement,  we  are  following  the  literature 
here) . 


and 


Accordingly,  ve  write 
U = f z1  (F*)"* 

a = h z (F*) 


(27) 


(28) 


where  f and  h are  dimensionless  functions,  hopefully  constants. 
We  must  have,  from  Fquation  27 


[ L j 1 [L4T~2]^  = L/T 
which  requires  j = ^ and  i = -1 . 
So 


U = f V/F*  /z 


From  Equation  28 

[L]  = [L]kfL4T‘2]£ 

which  requires  i = 0 and  k = 1 . 
So 

a = hz 


(29) 


(30) 


(31) 


(32) 
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and  furthermore,  from  Equation  19,  h must  be  a constant  equal  to  a. 

a = az  . (33) 

These  results  are  more  general  than  indicated.  Generalization 

to  a non-spherical  or  non-Hill  vortex  only  changes  the  numerical  value  of 

2 

a and  the  factor  - in  Equation  26. 

To  summarize,  we  have  found  for  any  time  similar  vortex  for  which 
assumptions  1,  2,  3,  and  4 hold,  the  radius  is  a linear  function  of  z and 
the  circulation,  which  is  proportional  to  Ua,  is  constant.  These  are  most 
important  findings. 

The  above  discussion  j ; in  no  sense  a proof  of  the  constancy  of 
circulation  or  the  linear  growth  of  vortices,  dimensional  analysis  can 
only  be  a guide  to  intuition.  Such  proof  as  exists  for  these  two  findings 
rests  on  laboratory  data.  However,  within  the  limits  of  error  of  labor- 
atory data  both  conclusions  are  verified,  although  one  could  wish  the 
laboratory  data  were  more  accurate. 

5.0  QUASI-PHYSICS  ON  VORTEX  SHAPE 

In  this  discussion  we  consider  the  relationship  between  energy 
and  impulse  in  an  isolated,  stable  vortex.  We  will  draw  some  inferences 
with  regard  to  the  likely  shape  of  an  atmospheric  vortex,  such  as  a rising 
fireball  or  thermal.  Our  inferences  cannot  be  firm  since  no  general  solu- 
tion exists  for  such  vortices,  thus  we  will  rely  heavily  on  intuition.  Our 
procedure  will  be  to  begin  the  discussion  with  the  Hill  spherical  vortex, 
the  best  understood  stable  vortex,  then  use  intuition  to  write  approximate 
formulae  for  other  .-hapes,  based  on  those  for  the  Hill  vortex. 
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We  will  base  our  discussion  upon  the  energy  of  motion  of  the 

vortex, 


dE 

dz 


= F 


(34) 


where  E is  energy,  F is  the  net  force  on  the  system,  and  z is  the 
distance  moved.  Here  we  have  tacitly  assumed  a time  similar  shape  in 
order  to  be  able  to  identify  single  characteristic  values  of  E,  F,  and 
z. 


Our  other  basic  parameter  is  the  impulse  uf  the  system, 


or 


(35) 


(36) 


where  P is  impulse  and  U is  the  characteristic  velocity  of  the  vortex, 
decal  1 that  the  impulse  of  a vortex  system  is  defined  to  be  the  momentum 
of  the  material  associated  with  vortex  motion  plus  an  amount  which  must 
be  radiated  to  infinity  whenever  the  momentum  is  changed. 


Hill  Spherical  Vortex 

For  a spherical  vortex  the  amount  of  momentum  lost  to  infinity 
is  j pVU  where  p is  the  exterior  density  and  V is  the  volume  of  the 
sphere.  Thus  the  impulse  of  a spherical  vortex  is  just 

P = 2ttp  a3U  , (37) 
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where  a is  the  radius  of  the  sphere  of  material  which  accompanies  the 
vortex  in  its  motion. 

The  energy  of  a Hill  spherical  vortex  is 

E = -y-  TTpa’u2  . (38) 

This  energy  is  composed  of  three  contributions,  (1)  the  energy  in 
the  exterior  flow  field, 

EX  = j (iTrpa3Lj2)  ’ (39) 

consistent  with  the  momentum  radiated  to  infinity,  (2)  the  energy  associated 
with  the  motion  of  the  center  of  mass  of  the  sphere, 

Ec»  = l(f’,pa3u2)  ' (A<» 

and  (3)  the  energy  associated  with  the  circulation  of  the  material  in 
the  spherical  vortex, 

Ep  = j 7Tpa3U2  . (41) 

Of  the  Equations  37  through  41,  only  38  and  41  are  peculiar  to 
the  Hill  vortex,  the  others  depend  only  on  the  spherical  shape  of  the 
vortex  (plus,  of  course,  the  assumptions  of  p*.o  , incompressibility, 

cl 

and  time  similarity). 

To  investigate  whether  an  expanding  vortex  with  flow  field 
instantaneously  identical  to  that  of  a Hill  vortex  is  an  appropriate 
approximation  for  vortex  motion  subject  to  an  external  force,  we  evalu- 
ate the  energy  requirements  for  a Hill  vortex.  From  Equations  34  and 
38, 
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(42) 


F 

But  we  know  that 


dz  7 \a  dz 


Ua  is  a constant  so 


2 dU  \ 
U dz  / 


dU  = _ 1_  da 
U dz  a dz 


and  Equation  42  becomes 


(43) 


dE  = 10H  pa3U2 


dz 


>/l  da\  _ E 
\a  dz  / a 


da 

dz 


(44) 


The  rate  of  delivery  of  energy  to  the  system  in  order  to  change 
the  momentum  is  obtained  by  substituting  Equation  37  into  Equation  36 


_ „ dP  _ n 3, ,2/3  da  1 dU\ 

F 1 u 37  - 2"Pa  u (a  d7  * u 37) 


(45) 


using  Equation  43  again  we  find 


U 4—  = 4Trpa3 
dz 


<£) 


(46) 


Comparison  of  Equation  46  to  44  shows  that  much  more  energy  is 
generated  in  changing  the  impulse,  by  whatever  force,  than  can  be  absorbed 
into  the  flow  field  of  a Hill  vortex. 

Remembering  that  we  have  radiated  j MU  of  momentum  away  to 
infinity  we  may  argue  that  all  but  j MU  of  the  energy  generated  in 
Equation  46  need  remain  in  the  system.  Thus  the  energy  which  must  be 
deposited  somewhere  in  the  flow  field,  E^,  may  be  as  small  as 


j Ttpa3U 


(47) 
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We  need  to  account  for  an  energy  excess  of 


dE  , 
d 

dz 


dE 

dz 


7Tpa3U2 


47  E da 
30  a dz 


(43) 


There  are  two  obvious  ways  out  of  this  dilemma.  The  first  is 
to  assume  the  excess  energy  goes  into  turbulence  and  the  second  is  to 
assume  a flatter  shape  for  the  vortex,  one  which  requires  more  energy. 
Both  of  these  are  plausible  and  indeed  have  data  to  support  them.  It 
is  known  that  buoyant  vortices  are  more  flat  on  top  than  a sphere  and 
in  some  cases  evidence  of  turbulence  is  observed. 


If  all  the  excess  energy  went  into  turbulence  it  would  require 
more  than  five  times  as  much  energy  in  turbulence  (plus  thermal)  than  in 
internal  circulation  energy  (47/9  to  be  exact,  from  Equations  41  and  48). 
It  is  most  difficult  to  imagine  that  a vortex  could  maintain  as  clean  and 
stable  a shape  as  is  generally  observed  in  the  laboratory  and  field  with 
five  times  as  much  random  energy  as  directed. 

Energy  Assumption 

On  the  basis  of  intuition,  guided  by  laboratory  and  field 
observation  we  will  assume  that  turbulent  energy  is  usually  negligible 
compared  to  circulation  energy  and  is  at  most  equal  to  circulation  energy. 

Consequences  for  Vortex  Shape 

The  energy  assumption  we  have  adopted  implies  a shape  consider- 
ably different  from  that  of  a Hill  spherical  vortex. 


No  convenient  equations  exist  by  which  one  can  accurately  study 
the  dependence  of  energy  on  vortex  shape,  at  least  none  such  are  known  to 
this  author. 
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Two  cases  are  available,  the  Hill  spherical  vortex,  with  a 
continuous  smooth  vorticity  distribution,  and  that  of  a ring  of  uniform 
vorticity  with  cross  section  small  compared 
to  its  distance  (see  figure) . Based  on  the 
relations  on  p 241  in  Lamb1  the  circulation 
energy  of  such  a vortex  can  be  many  times 
larger  than  the  energy  of  translation.  We 
accept  this  as  an  indication  that  a flat- 
tened vortex  will  result  in  increased  cir- 
culation energy  relative  to  translations^ 
energy  and  sufficient  flattening  will  in  fact  increase  the  circulation 
energy  to  any  degree  we  require. 

Our  next  task  is  to  quantitatively  estimate  the  increase  of 
energy  with  flattening.  Since  we  have  no  other  convenient  starting  point 

we  will  estimate  the  likely  effects 
beginning  with  Hill's  spherical  vor- 
tex of  radius  a.  The  solid  curve 
represents  a flattened  spheroidal 
vortex  with  horizontal  radius  a' 
and  vertical  radius  b.  The  volume 
is  to  be  the  same  for  both  cases  , 

y TTba'2  = y TTa3  (49) 

and  the  rise  velocity  is  to  be  fixed 
'J'  = U (50) 

as  is  the  mass. 
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For  convenience  define 


f = b/a  , 


(51) 


then 


a'  = a/Vf-  . 


(52) 


We  expect  the  impulse  of  the  new  vortex  to  be  about  equal  to 
that  of  the  Hill  vortex, 

P'  = momentum  of  moving  material  + radiated  momentum 
proportional  to  volume 

= MU  + j MU  = 2irpa3U  = P (53) 

since  ba^2  = a3  , from  Equation  49. 

The  energy  associated  with  the  center  of  mass  motion  is  also 

unchanged 


ecm  " f "oa>u2 


(54) 


The  energy  associated  with  circulation  should  be  increased  due 
to  the  necessity  for  the  circulating  material  to  go  farther  out,  in  the 
ratio  a'/a,  and  to  have  less  time  to  get  there,  in  the  ratio  b/a.  This 
time  factor  is  required  because  we  must  have  continuity  of  tangential 
velocity  at  the  interface.  Thus  the  internal  velocity  of  circulation 
should  be  increased  by  a factor  of  roughly 


a^a 
a b 


a 7b 


(55) 


and  the  corresponding  energy  increased  by  the  square  of  this  factor, 
from  Equation  41  we  write 

Ep  = j Tipa3U2(a7b)2  = j ^pa3U2/f3  . (56) 

The  energy  external  to  the  vortex  should  also  be  increased 
since  material  near  the  vortex  surface  must  also  flow  around  the  vortex 
faster  in  less  time.  Almost  all  of  the  energy  external  to  the  Hill  vor- 
tex resides  in  material  very  close  to  the  edge  because  the  velocity  falls 
off  as  the  third  power  of  the  distance  (thus  energy  as  the  sixth  power) . 
Perhaps  as  good  a way  as  any  to  estimate  the  new  energy  is  to  assign 
the  same  relative  increase  in  velocity  as  that  for  the  internal  motion 
to  the  same  amount  of  material  as  was  the  case  for  a Hill  vortex.  This 
process  yields,  from  Equation  39 

E ' = Ex(a7b)2  = | pa3U2/f3  (57) 


The  energy  of  the  flattened  vortex  is 


E " = E".  + Er"  + E'  = 
CM  f x 


10 

7 


Trpa  3U  2 


[l5  * 15f  > ] 


= E (7  + 8/ f 3 ) / 1 5 


(58) 


Comparison  of  Equation  58  to  Equation  47  shows  that  we  require 


(7  + 8/ f 3 ) / 1 5 = 77/30 


(59) 


if  no  appreciable  amount  of  energy  causes  turbulence.  For  this  esse 

f ~ 0.633  , or  b/a'  = 0.506  (60) 

Note  that  the  ratio  b/a'  will  be  the  observed  ratio  of  vertical 
to  horizontal  radius. 
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If  we  take  the  other  extreme  of  our  energy  assumption  we  can 
assume  an  amount  of  energy  equal  to  E'  is  in  turbulence,  then  Equa- 
tion 58  becomes 

E'  = E'm  + 2E'  + E'  = y-  irpa3U2  (14  + 25/f3)/30  (61) 

We  then  require,  from  Equation  47,  that 

(14  + 25/f 3)/30  = 77/30  (62) 

which  has  the  solution 

f = 0.737  , b/a'  = 0.642  (63) 

Turner3  on  page  188  reports  Scorer  as  having  found  a value 
somewhat  above  our  range,  b/a'  = 0.716  from  laboratory  experiments 
dropping  salt  water  "thermals"  through  pure  water.  In  Turner's  book 
a shadowgraph  is  illustrated  which,  by  this  author's  measurement,  yields 

0.49  < b/a'  < 0.56  , (64) 

in  good  agreement  with  our  prediction  (Turner's  shadowgraph  shows  clear 
evidence  of  turbulence  , but  not  extreme  turbulence).  Similar  measurements 
on  Figures  10  and  15  of  Reference  5 yield  b/a'-  .65. 


Until  more  definitive  data  is  uncovered  (it  almost  certainly 
exists),  we  will  conclude  that  vortices  under  an  external  force  are  not 
very  spherical,  but  can  be  expected  to  be  flattened  spheroids  with  ratio 
of  vertical  to  horizontal  radius  of  about  0.6. 
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6.0 


VORTEX  IN  UNIFORM  MEDIUM 


In  this  section  we  derive  the  useful  results  of  Section  4 for 
a flattened  spheroidal  vortex. 

The  Equation  17  is  still  valid, 

dV  = — dM  , (65) 

Ka 

and  we  will  continue  to  write  the  mass  entrainment  equation  as 
dM  2 

HLi  = 47Ta  p a , (66) 

dz  a v 1 

where  a will  henceforth  be  understood  to  mean  the  horizontal  radius. 
The  vertical  radius  will  be  written 


b = f a 
v 


(67) 


If  we  carry  out  the  differentiation  of  volume  with  respect  to 
altitude  in  Equation  65  and  use  Equations  66  and  67  we  find 


. r 2 da  1.0 
4tt  f a^  -r-  = — 4tt a^p  a 
v dz  p ya 

' a 


(68) 


or 


da  .. 

A~  ~ a/f 

dz  v 


(69) 


This  is  the  only  effect  of  aspherical  shape  on  the  results  of 
Section  4.  In  particular,  the  buoyancy  is  still  constant,  as  is  the 
quantity* 

X"  E Ua  (70) 

* We  use  X'  in  preference  to  circulation.  They  differ  only  by  a 
numerical  constant. 
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If  we  rewrite  Equation  25  the  factor  f cancels  and  the 
momentum  equation  becomes 

ft  (a3U)  = § g (pa  - p)  a3/ pi)  (71) 

where  we  have  written  — as  U -r- 

dt  dz 


Remembering  that  the  quantity  X'  is  constant  we  can  write 
Equation  71  as 


_d_ 

dz 


2g(Pa  - P)  a3 
p3)Pu 


(72) 


Carrying  out  the  indicated  differentiation  in  Equation  72,  sub- 
stituting Equation  69  and  rearranging  gives 


g(pa  - p)  a3fy 

a _ 

3p  X"2 


(73) 


We  note  that  the  numerator  is  proportional  to  the  buoyant 
force  which  is  constant,  p^p  , and  X'  is  constant.  Our  result  is 
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therefore  consistent  with  our  earlier  conclusion  that  a is  a constant 
and  contains  the  surprising  information  that  the  entrainment  varies 
inversely  as  the  square  of  the  circulation  (or  X'),  that  is,  the  more 
circulation  present,  the  less  entrainment. 


This  dependency  is  said  to  be  verified  in  laboratory  experi- 
ments where  the  assumptions  of  incompressibility  and  p - p are  very 
good,3  and  the  circulation  can  be  varied.  We  have  personally  witnessed 
demonstrations  which  qualitatively  verify  Equation  73.  In  addition,  the 
above  treatment  is  good  within  experimental  error  for  atmospheric 
thermals3  which  in  nature  are  weakly  buoyant  (p  ^ p ) and  rise  a dis- 
tance  small  compared  to  a scale  height. 
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Encouraged  by  the  above,  we  will  continue  to  extend  the  treat- 
ment toward  the  case  of  interest. 

7.0  EXTENSION  TO  ADIABATIC  LAPSE  CASE 


Our  first  step  will  be  to  make  a "plausible"  extension  to  a 
real  atmosphere  which  has  an  adiabatic  lapse  rate.  We  are  relaxing 
the  assumption  of  incompressible  flow  but  keeping  the  approximation 
p = p and  arranging  a situation  where  buoyancy  is  constant. 

3. 

For  convenience,  define  the  density  factor  to  be  given  by 


A E pa  - p 


Then  Equation  71  becomes 


A roa 3 m = 2gajA 

dz  tp  J 3U 


The  entrainment  equation  becomes 


-f-(pa3)  = 3p  aza/ f 
dzv  J a v 


(74) 


(75) 


(76) 


The  volume  equation  requires  two  terms 


+ 

adiabatic 


entrainment 


(77) 


Under  any  of  the  three  assumptions,  (i)  incompressibility,* 
(ii)  no  small  scale  mixing,  or  (iii)  pressure  equilibrium  where  interior 
and  exterior  gas  constants  pf)  are  equal,  the  entrainment  term  is  given 


* No  longer  true  for  this  case. 
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by 


i“3 


= (T  i (pa!)  = 3aVfv  • 

a 


entrai ament 


(78) 


where  we  have  used  Equation  76. 


The  adiabatic  term  is  given  by 


a3  dp 

, . , vp  dz 

adiabatic 


(79) 


where  p is  ambient  pressure,  assumed  equal  to  vortex  mean  pressure. 


The  volume  equation  is  thus 


_d_ 

dz 


JL3  + 
YP  dz 


3a2a/f 

v 


(80) 


The  only  equation  which  has  a different  appearance  from  the 
incompressible  case  is  Equation  80.  We  will  first  show  that  the  buoy- 
ancy term  in  Equation  75  is  constant 


_d_ 

dz 


(a ’A)  . 


' Jl  (pa>) 


(81) 


Substituting  Equations  76  and  80  we  find 


_L  dp 

yp  dz 


(82) 
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But  the  condition  for  an  adiabatic  lapse  rate  is  that  the  "buoyancy 

gradient" 


B 0 

Pa  dz  Yp  dz 


So  we  find  indeed  that  the  buoyancy  is  constant 


± (.*4)  - 0 


We  can  use  this  fact  to  write  Equation  75  as 


£ (Pa’u)  - 


2«V. 

3U 


where  the  subscript  o indicates  initial  value. 


(83) 


(84) 


(85) 


At  this  point  we  have  used  up  Equation  80,  having  found  it  to 
be  equivalent  to  Equation  76, and  find  ourselves  with  two  equations  in- 
volving the  quantity  pa3,  which  look  entirely  similar  to  the  incom- 
pressible case.  The  dimensional  arguments,  unconvincing  at  best,  which 
led  us  in  the  analysis  of  the  incompressible  case  are  less  useful  now 
since  we  have  an  additional  dimensioned  variable,  p.  We  will  therefore 
adopt  two  methods  of  reasoning  by  analogy,  one  essentially  mathematical 
and  one  more  physical.  The  results  will  be  identical. 


The  Mathematical  Guess 


We  will  assume  that  the  quantity  to  be  conserved  is  the 
product  of  U and  the  cube  root  of  its  coefficient  in  Equation  85, 

X = p1'3  X"  = pl/3  a U * const.  (86) 
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There  is  no  justification  short  of  experimental  comparison 
for  this  assumption,  as  was  the  case  for  X'  = const,  in  incompressible 
flow.  We  defend  it  here  only  on  the  basis  (1)  it  seems  as  reasonable 
as  most  alternatives,  (2)  the  consequences  of  the  assumption  are  simple 
and  seem  more  reasonable  than  most  other  alternatives  we  have  tried. 


Rewrite  Equation  85,  taking  into  account  the  constancy  of 


X, 


d . 3.2/3  “gAoao 

dl  (pa  } 


3UX 


(87) 


Carry  out  the  indicated  differentiation  to  find 


2g4«a»  fV 

2p  a “ ' -^Ux 


(88) 


where  we  have  used  Equation  76.  Solving  for  a yields 

gA.a]  f 
& 0 0 v 
a = 

3p1/3  X2 


(89) 


We  find  a to  be  no  longer  constant  but  to  vary  as  the  in- 
verse cube  root  of  the  density. 


Substitution  of  Equation  89  into  Equation  76  now  allows 

solution, 


(90) 


(91) 
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Substitution  of  Equation  91  into  Equation  85  yields 


u = (92) 

The  above  treatment  seems  to  yield  reasonably  satisfactory  re- 
sults, intuitively,  but  one  is  hesitant  to  base  an  argument  on  the  assump- 
tion that  the  third  moment  of  the  velocity  distribution  is  constant.  That 
is, . 


X3  = pa3U3«  MU3  = const.  (93) 

While  this  is  fascinating  and  may  be  true,  one  would  prefer  to 
derive  it  rather  than  assume  it. 


The  Physical  Guess 


As  an  alternative  way  of  thinking  about  the  problem  we  consider 
the  equation  which  introduces  the  entrainment  parameter,  a,  and  solve  it 
for  a in  terms  of  dynamical  quantities.  Solving  Equation  76  we  find 


a 


a dM 
M dz 


(94) 


where  M is  the  mass  associated  with  the  motion.  But  in  the  incompres- 
sible case  we  found  a could  be  expressed  as,  from  Equation  73 


a 


gA „ a „ f 
6 oo  v 

3pa2U2 


a dE 
E dz 


(95) 


where  we  have  taken  M~pa3  and  &~pa3U2. 


We  imagine  that  the  forces  operating  on  the 
relation  between  mass  entrainment  defined  by  Equation 


fluid  require  a 
94  and  energy  buildup 
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defined  by  Equation  95  to  balance.  Then  a may  always  be  given  by  an 
expression  like  Equation  95,  whether  flow  is  incompressible  or  not  or 
indeed  even  in  cases  where  buoyancy  is  not  constant.  Thus  we  will  assume 


a = f 

V 3pU2 


(96) 


Differentiation  of  the  left  hand  side  of  Equation  85  and  use  of 
Equations  76  and  96  along  with  84  yields 


d , 3...  8V.  3 dU  2SV 

57  tpa  u)  = — pa  " 


dz  3U 


(97) 


or 


-r  * * 3 

& 0 0 


(98) 


Substitution  of  Equation  96  into  Equation  76,  multiplication  by 
U2  and  making  use  of  Equation  84  yields 


U2  4-  pa3  = gA  a3 

uZ  0 0 


(99) 


Addition  of  Equations  98  and  99  shows 


a?  (pa3u2)  = gVo/3 


(100) 


SO 


pa3U2  = p a3  U2 
ooo 


+ \ 8^oao  z 


(101) 


Solve  Equation  101  for  U2,  substitute  into  Equation  99  and 
rearrange  to  find 
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1 


pa 


r£(pa3)  = 


poUo  l l-i 

gA0  + 3 ZJ 


(102) 


which  has  the  solution 


a = a 


/po\1/3/ 

At)  1 + 


§Ac 


3poU' 


\ 

z; 


(103) 


Equation  103,  derived  from  physical  intuition,  is  idential  to 
Equation  91  which  was  based  on  mathematical  intuition.  It  follows  that 
Equations  92  and  86  are  also  true  under  this  treatment  . 


8.0  EXTENSION  TO  REALISTIC  CASE 

In  the  general  case  of  a real  fireball  vortex  in  a real  atmos- 
phere processes  occur  which  are  certain  to  alter  the  circulation  and  X. 

The  most  straightforward  approach  is  to  attempt  to  list  all  such  processes, 
estimate  the  effect  on  X,  then  follow  through  the  consequences  without 
altering  any  other  assumption.  This  will  be  our  approach.  We  will  find 
additional  terms  in  the  equations  for  X,  a,  buoyancy,  etc.  These  will 
form  the  set  of  equations  to  use  in  comparison  to  data. 


8.1  Dry-Air  Lapse  Rate  Effect  on  X 


At  a particular  moment  in  the  history  of  the  vortex  it  will  be 
entraining  air  which  originally  resided  at  an  altitude  not  equal  to  the 
vortex  altitude.  It  will  turn  out  to  have  originated  mostly  from  below 
the  vortex  altitude.  The  effect  on  circulation  can  be  calculated  via 
Bjerknes'  Theorem: 


(104) 
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where  the  partial  derivative  indicates  we  are  concerned  only  with  the 
change  in  circulation  due  to  lapse  rate. 

If  we  assume  pressure  equilibrium,  carry  the  integral  up  the 
axis  from  to  + °°,  and  complete  the  circuit  in  a region  where  motion 
is  negligible,  we  find 


T 


where  the  limits  include  the  region  of  disturbed  density  (i.e.,p  t p ). 
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The  value  of  p can  be  determined  by  using  the  fact  that  each  air  parcel 
had  ambient  values  pfl,  p , and  has  expanded  adiabatical ly  to  pressure 
equilibrium,  p,  at  the  current  vortex  position. 


P * P0(P/P0)1/Y 


- 52 


dz, 


i/Y 


Expanding  to  first  order 


(106) 


1_ 

P 


1_ 

P, 


B 


(107) 


where,  as  before, 


I Je. . J_  3s. 

p cte  ' yp  9z  ’ 


(108) 


is  the  buoyancy  gradient  of  the  ambient  atmosphere  am.  5 z'  is  the  mean 
distance  below  the  fireball  from  which  the  material  along  the  axis  origi- 
nated. The  rate  of  change  of  the  circulation  becomes 
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ar 

dz 


L 


T 

i ar  _ vp  f 52"Bd£ 
u at  l “ " u J p 


(109) 


B 

It  is  difficult  to  know  what  limits  to  assign  to  the  integral 
for  a realistic  case.  The  upper  limit  is  clearly  near  and  a bit  beyond 
the  top  of  the  sphere,  call  it  distance,  a,  above  the  center,  the  lower 
limit  might  be  as  low  as  z-  5z  since  all  of  the  material  above  the  point 
will  eventually  be  entrained  into  the  vortex  and  the  material  below  will 
not.  More  sophisticated  reasoning  might  lead  to  a better  answer  but  we 


will  accept  these  limits  for  the  present 


z+a 


z-6z 


(HO) 


For  simplicity,  we  will  assume  5 z«H  and  that  the  buoyancy 
gradient  is  constant  over  the  interval,  thus 


z+a 


z+a 


dr 
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= . M f 

L Up  J 
= -fi  = 


5 z'd£ 


-U 


5 z 'di 


z-<$: 


^ gB<5z2 
2U 


We  have  also  assumed  5 z'«i. 


(Ill) 


In  a later  report  we  will  address  the  question  of  the  value  of  5 z, 
which  depends  on  the  entrainment  parameter,  a,  and  the  shape  of  the  vortex. 


Converting  Equation  111  into  an  equation  for  X yields 


jPC  _ gp2/3a6z2 
3z  L " 2CX 


(112) 
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Our  justification  for  carrying  p1/3  through  the  derivative 
is  that  we  have  already  assumed  adjustments  in  the  flow  field  to  occur 
which  cancel  density  variations  for  B=0,  in  our  discussion  of  the 
adiabatic  lapse  rate  case.  The  constant  C relates  T to  Ua,  i.e.,  r = CUa. 

8.2  Effect  of  Internally  Generated  Buoyancy 

When  dissociated  fireball  atoms  combine  to  form  molecules  heat 
energy  is  released,  which  increases  the  buoyancy  of  the  region.  Water 
vapor  previously  entrained  from  the  ground,  lagoon,  or  along  with  ambient 
air  can  also  condense,  releasing  considerabi  quantities  of  heat  energy. 

These  cases  appear  much  more  difficult  to  handle  than  that  just 
discussed.  The  conceptual  problem  arises  because  release  of  dissociation 
or  condensation  energy  may  take  place  partly  or  wholly  in  regions  remote 
from  the  axis.  Thus  there  is  no  obvious  effect  upon  either  Vp  or  p 
along  the  axis  and  therefore  no  obvious  effect  upon  the  circulation.  On 
the  other  hand,  one's  intuition*  is  repelled  by  the  concept  that  a sudden 
release  of  buoyancy  at  t=0  will  cause  the  generation  of  a great  deal  of 
circulation  whereas  a release  of  buoyancy  later  will  not  affect  the  cir- 
culation, provided  the  region  of  release  is  a centimeter  or  two  from  the 
axis,  independent  of  the  amount  of  buoyancy  released. 

The  paradox  can  be  resolved  qualitatively  be  remembering  that 
the  vortex  moves  as  a unit,  held  together  by  dynamic  forces  which  we  may 
not  understand  at  all  well,  but  which  nevertheless  succeed  in  holding 
real  vortices  together  in  the  real  world.  With  this  point  in  mind  it 
seems  clear  that  creation  of  a buoyant  element  anywhere  in  the  vortex 


*At  least  this  author's  intuition. 


37 


causes  creation  of  pressure  gradients  which  do  not  destroy  the  integrity 
of  the  vortex  (experimental  fact)  but  which  must  cause  the  vortex  to 
accelerate  upwards  precisely  according  to  the  amount  of  buoyancy  released, 
averaged  over  the  entire  vortex.  If  we  could  calculate  these  gradients 
we  would  be  able  to  apply  Bjerknes'  theorem  to  find  the  rate  of  increase 
ot  circulation.  Since  we  have  not  yet  succeeded  in  solving  the  problem 
we  are  forced  to  use  a simple  argument  which  we  hope  yields  the  correct 
dependency  but  which  is  likely  to  be  wrong  by  a numerical  factor. 

We  will  assume  that  only  the  change  in  average  density  due  to 
additional  buoyancy  is  effective  in  circulation  generatior  (we  could 
reach  the  same  answer  using  an  equivalent  argument  based  on  Vp) . The 

change  in  density  caused  by  release  of  energy  6 E followed  by  adiabatic 
expansion  is 


Vp  = -1  5E 
P YE 


(113) 


The  residual  in  Bjerknes'  theorm  is,  since  we  assume  all  other 
factors  must  cancel, 

f = - g f fiz  (114) 

B B 

Now  the  question  arises  as  to  how  long  the  circulation  genera- 
tion is  operative.  We  know  that  when  the  vortex  was  originally  generated 
from  a buoyant  bubble  at  rest,  circulation  increased  during  roughly  the 
period  required  for  the  material  to  move  up  the  axis,  i.e.,  during  the 
time  required  for  the  vortex  to  turn  through  a radian  or  two.  Since 
our  i.ieihanism  is  not  essentially  different  we  can  expect  the  appropriate 
time  to  be  given  by 


T 


(lib) 
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Multiplication  of  Equation  114  by  Equation  115  yields 


5 r 


2a2g  5E 
yU  E 


(116) 


or 


Pv  p1/35  T ^ 2a3p2/3  g SE 
dA  " C CyX  E 

where  C is  the  factor  which  relates  T to  Ua  . 


(117) 


Replacing  the  5 's  with  derivatives  with  respect  to  z yields 


8.3 


ax 

3z 


2gp2/3a3  1_  dE 
H‘  CyX  E dz 


(118) 


Effect  of  Wind  Shear  on  X 


On  the  average,  the  vortex  should  be  transported  with  the  mean 
horizontal  (and  vertical)  velocity  of  the  local  wind.  However,  vertical 
gradients  in  the  horizontal  wind  velocity  can  be  expected  to  distort  the 
vortex  and  destroy  circulation. 


The  cross  hatched  area  in  the 
sketch  illustrates  the  effect  in  question. 
Vorticity  transported  across  the  center- 
line  will  have  the  wrong  sign  on  the 
new  side  and  will  cancel  an  equal  amount 
of  vorticity.  Thus  if  the  vorticity 
were  uniformly  distributed  we  would  have 


5r  = 


>/2\  i_e 

\ 7T  / 77 


(119) 


39 


where  the  factor  of  2 arises  from  the  loss  of  vorticity  on  one  side  plus 

2 

cancellation  on  the  other  side  and  the  factor  — arises  from  an  azimuthal 
average  of  the  cross  product  of  the  vorticity  vectors  (<cos2<J>>) . The 
angle  5 0 is  simply  given  by 


dO  _ 1 dv^ 
dz  U dz 


(120) 


The  switch  from  T to  X is  in  this  case  direct,  substituting 
A for  I'  in  liquation  119  and  X for  U in  Equation  120  yields 


9X 

3z 


(121) 


We  should  not  be  surprised  if  Equation  121  yields  somewhat  too 
large  an  effect  since  the  vorticity  is  likely  to  be  concentrated  away 
from  the  axis. 


8.4  Equation  for  X 


If  we  define 


e ? P1/3  a 


(122) 


After  the  vortex  is  formed  we  anticipate  an  equation  for  X 
which  has  the  form 


dX 

dz 


- 
Sz  L 


ax 

dz 


3X 

az 


= I dE  \ £C  IZj 

CX  \2a  y E dz  / 7i^v  \ 3Z 


M2 


8v 


(123) 
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Unfortunately,  at  the  present  state  of  knowledge,  each  of  the 
terms  in  Equation  123  may  be  expected  to  be  in  error  by  a multiplicative 
factor  of  order  unity  (at  best). 

8.5  Equations  for  a,  £ 


The  momentum  equation  (75)  becomes 


A 

dz 


(C2x) 


2«a3A 

3X 


(124) 


and  the  entrainment  equation  (76)  becomes 


p 1/3  a/  fy 


(125) 


If  we  carry  out  the  indicated  differentiation  in  Equation  124  and  sub- 
stitute Equation  126  we  find 


f ga3A  f £ lv 

v6  _y_L_  dx 

3p1/3  X2  2p*'3X  dz 


(126) 


The  first  term  is  the  same  as  we  found  for  incompressible 
flow  and  assumed  for  compressible  flow  in  an  adiabatic  lapse  rate, 
except  that  now  a3A  may  vary.  The  second  term  is  required  by  momen- 
tum balance  provided  our  assumption  for  the  adiabatic  lapse  rate  case  is 
correct  , and  time  similarity  is  still  a good  assumption. 


Substitution  of  Equation  126  into  125  yields  the  equation 

for  £. 


A f = . £ £X  (127) 

dz  " 3X2  2 X dz 
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8.6 


Equation  for  Mass  Discrepancy 


Except  for  the  term  in  the  mass  discrepancy,  Aa3,  the 
Equations  123  and  127  form  the  complete  set  necessary  to  solve  a problem. 
We  now  write  the  equation  for  the  mass  discrepancy  term, 


Jj(a’A)  “ E(a  Va!p) 


dPl 

dz 


da  ,3 

p -j 3a  p a/f 

a dz  a v 


(123) 


the  last  term  is  obtained  from  the  entrainment  Equation  76. 

The  derivative  of  volume,  or  a3,  has  three  contributions, 

(129) 


where 


d 3 9 3 

j — aJ  = -r-  aJ 
dz  9z 


+ — a3 

+ — a3 

....  9z 
adiabatic 

, . 9z 

heat 

entrainment 


9 3 

-s-  aJ 
9z 


_9 

9z 


_9 

9z 


1 adiabatic 

1 a3  dE 

cl 

ihoatYE  dz 

a3 

a3  dp 
yp  dz 


= 3a  ^ a/f 


(130) 


(131) 


(132) 


entrainment 


Substitution  of  Equations  130,  131,  and  132  into  128  yields 

3F  (a!«  - «’(B  * £ af)  • (133) 

This  completes  the  set  of  equations  necessary  to  solve  the  vor- 
tex motion,  given  initial  conditions  on  X,  £,  and  a3A. 
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8.7 


Removal  of  the  Assumption,  p > pa 

a 

The  careful  reader  will  note  that  in  the  equations  defining  X 
(86)  and  £ (122)  we  used  the  symbol  p with  the  understanding  that  it 
could  be  regarded  either  as  mean  vortex  density  or  as  ambient  density,  the 
two  being  considered  near  enough  to  equality. 

If  one  actually  applies  this  theory  numerically  one  must  decide 
which  density  to  use  to  initialize  X and  £,  and  do  it  at  zero  time 
when  the  densities  are  very  far  from  equal.  In  practice  we  have  used  p 

cl 

for  both  of  these  quantities  but  one  wonders  whether  the  early  stages  of 
the  calculation,  when  p f p cause  a significant  effect  on  the  final 
answer. * 


The  most  straightforward  method  of  answering  this  question  i.:> 
to  remove  the  approximation  p ~ p#  . We  proceed  to  do  this.  So  far,  in 
practice  it  has  not  yielded  significantly  different  results  from  previous 
equations  embodying  the  approximation.  We  preseat  it  for  completeness. 

Define 

£ = p1/3  a , X 5 p1/3  all  , 8 = a3A  = (p  -p)a3  . (134) 

cl  3 3 

The  entrainment  equation  (76)  becomes 

^(£3-B)  = 3P"3  £2  a/f v , (135) 


* Beyond  a timing  error  comparable  to  torus  time  and  an  altitude 
error  comparable  to  the  initial  value  of  a. 
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or 


d£ 

dz 


a/*v 


i da 

Ic2"  dz 


and  the  impulse  equation  becomes* 


ElO*  Pa/2>a3uI 


ga3A 

l! 


(136) 


(137) 


(138) 


Carry  out  the  indicated  differentiation  in  Equation  138  and 
make  use  of  Equation  135  to  solve  for  et 


a 


(139) 


When  a/C3  = A/p  is  small,  the  first  two  terms  on  the  right 
of  Equation  139  yield  the  previous  result,  Equation  126.  It  is  not  ob- 
vious that  the  third  term  is  negligible  but  estimates  for  atmospheric 
cases  indicate  it  to  lie  in  the  range  of  10  1 to  10  compared  to  the 
first  term. 


The  equation  for  a (139)  and  the  equation  for  £ (136)  are 

the  necessary  modifications  to  include  the  case  p / p^.  As  was  previously 
pointed  out,  we  have  not  found  a case  where  this  change  caused  a significant 
numerical  effect  in  any  observable  quantity. 


* See  discussion  in  Section  5 for  the  factoring  of  impulse 
into  interior  and  exterior  terms. 
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9.0 


INITIAL  CONDITIONS 


We  must  invent  some  way  to  start  a calculation.  The  simplest 
method  is  to  derive  values  of  a,  X,  £,  and  6 from  scaling  laws  or  simple 
theory  and  set  these  values  at  burst  time  and  altitude.  Thus  the  early 
motion,  presumably  up  until  about  torus  time,  is  treated  as  instantaneous. 
Clearly  this  is  inadequate  if  one  is  interested  in  early  times,  and 
straightforward  to  correct.  However,  if  one  is  interested  in  intermediate 
or  late  motion  it  introduces  an  error  in  altitude  of  at  most  the  order  of 
the  initial  radius  and  in  time  of  at  most  the  order  of  torus  time,  both 
small  perturbations. 

9.1  Free  Air  Shock  Circulation 

The  circulation  generated  by  reaction  forces  in  giving  off  a 
shock  by  a free  air  burst  in  an  exponential  atmosphere  has  been  estimated.6 

'Hie  formula  is 


r = 4.2  X 105  (sec-1)  aj/ll  , (140) 

where  aQ  is  pressure  equilibrium  fireball  radius  as  usual  and  H is 
atmospheric  e-folding  distance.  In  terms  of  X,  if  we  assume  a factor 
of  5 scaling  from  circulation  to  Ua,  we  find 

X'  = 8.4  x 104O1/3  a?/H  (141) 

r a o v ' 

We  have  placed  primes  to  indicate  that  these  values  of  F and 
X are  not  valid  for  all  values  of  afl. 
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(U)  The  formula  (141)  results  from  a first  order  expansion  treatment. 

The  largest  expansion  parameter  reaches  a value  of  0.1  when  ao/H  = 3.9  x 
lO-2  (or  aQ^  250  m) . For  larger  pressure  equilibrium  radii  the  formula 
becomes  unreliable.  We  need  to  invent  some  way  to  extend  the  formula. 


A simple  method  to  try  is  based  upon  the  fact  that  Equation  141 
will  yield  too  large  a value  for  mean  rise  velocity  when  aQ~H,  by  a 
factor  of  3 or  4.  We  might  try  formula  (141)  until  a0  is,  say  twice  as 
big  as  the  value  quoted  above,  then  simply  hold  rise  velocity  constant 
for  the  next  factor  of  3 or  4 in  aQ.  That  is,  we  take 


X 


r 


8.4  x I04p1/3 
a 


8.4  x 10 4p1/3 
a 


V./H 


a <a 

o 1 


a <a  £ 4a , 
10  1 


(142) 


where 


5 x 


104  cm  or  4 km. 


This  argument  seems  fairly  reasonable  and  turns  out  to  give 
fairly  good  agreement  with  data  up  to  the  limit  4a j but  should  be 
verified  (improved). 
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9.2 


Reflected  Shock  Circulation 


No  very  convincing  calculation  of  the  effect  on  circulation  of 
a shock  traversal  of  a vortex,  or  even  a spherical  fireball,  exists.  We 
will  develop  a formula  which  might  be  worth  trying  but  which  should  be 
viewed  with  great  skepticism. 


Brian  Murphy  has  developed  a formula  for  the  displacement 
velocity,  V^,  of  the  edge  of  a spherical  pressure  equilibrium  bubble 
struck  by  a plane  shock.  The  formula,  in  our  units,  is 


V 


D 


C 

a 


A 

Pa+P 


f (x) 


(143) 


where 


x = i.135  zd  (p/nT)1/3  (144) 

is  ambient  sound  speed,  zd(cm)  is  distance  from  the  center  of  the 
image  explosion  (or  real  explosion  in  a multiburst  situation),  p (dynes) 
is  ambient  pressure,  (erg)  is  the  total  yield  of  the  image  explosion, 
and 


.0814  x~3  + ■ 2086x~  3/2 
[1  + . 096x~ 3 + . 2504x"3/2]‘/2 


(145) 


To  the  degree  of  approximation  possible  we  can  assume  f^Vp  is 
given  to  all  of  the  material  inside  the  fireball,  where  f ~1.  The  incre- 
ment of  momentum  delivered  to  the  fireball  material  is  then 


6P 


r 


4 3 _ , 

3 *pa  frVD 


(146) 
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By  the  time  the  vortex  flow  field  has  readjusted  outside  as  well 
as  inside,  this  increment  in  momentum  will  result  in  an  increment  in  vortex 
rise  velocity  given  by  momentum  conservation, 

j 7T (p  + pa/2)  a35Ur  = 5Pr  , (147) 


or 


6U  = 


f V = f C 


(P3-P)P 


G p+p  2 r D r a , 

M Ma  p +p 

a 


)Hk 


f(x)  . 


and  therefore  the  reflected  shock  contribution  to  X is 


(148) 


(P„-P)P 

5X  = f p1/3  C a - — f (x) 

G ra  a /pm/plp; 


(p-p)(k  ;j 


(149) 


The  formula  (149)  has  the  satisfactory  property  of  going  to  zero 
as  A ■+  0 but  the  unsatisfactory  property  of  going  to  zero  as  p/p-+-0.  In 
fact  5X  is  a maximum  with  respect  to  p/p  at  p/p  ~ 0.4. 

3 3 

Perhaps  this  is  reasonable  but  it  is  certainly  not  obvious. 


Another  unsatisfactory  aspect  of  Equat  n 149  is  that  it  diverges 
for  Zj  = 0.  The  circulation  of  a ground  burst  may  indeed  differ  from  that 
of  a free  air  burst  of  the  same  yield  but  it  is  not  at  all  obvious  diether 
it  is  greater  or  less.  That  is,  the  mass  and  buoyancy  are  about  the  same, 
a strong  reaction  force  due  to  the  ground  should  tend  to  add  circulation, 
but  the  development  of  the  flow  field  may  be  severely  hampered  by  the  pres- 
ence of  an  impenetrable  boundary. 


For  the  present  we  will  suggest  that  f(x)  should  be  held  to  a 
finite  value  or  reduced  to  zero  as  burst  height  becomes  less  than  a. 
Clearly,  this  is  a point  which  deserves  some  attention. 
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Digression  to  Multibursts 


Here  we  digress  to  a more  general  consideration.  The  ground 

effect  treatment  could  form  the  basis  for  a multiburst  shock  effect  if 

we  generalize  our  theory  slightly.  First  calculate  the  value  of  the 

shock  induced  circulation,  X , according  to  Equation  149  but  with  z. 

s d 

in  Equation  144  replaced  by  r,  where  r is  the  distance  between  the 
center  of  the  vortex  and  the  new  explosion.  Let  X be  a vector.  The 

“►  -f 

contribution  Xg  is  a vector  pointing  along  the  vector  r.  Add  Xg  to 
X vectorally.  Equation  123  becomes  the  equation  for  the  vertical  compo- 
nent of  X.  Whenever  X enters  an  equation  for  a,  £ , or  B,  use  |x|. 
The  equation  for  the  horizontal  component  of  X is  X = 0. 

We  have  noc  tried  this,  but  presumably  comparisons  could  be 
made  to  Reference  8 . 

9.3  Buoyancy  Circulation 

Calculations  have  been  carried  out  for  an  initially  spherical 
bubble  rising  under  buoyant  forces  for  the  extreme  case  of  zero  internal 
density,  pressure  equilibrium,  incompressible,  inviscid  flow.6  These 
calculations  yield 

rB"  4-5  a o V^"o  ’ (150) 


and 


U = 0.8 

B 


(151) 
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The  method  of  calculation  does  not  converge  all  the  way  to  torus 
time.  The  above  numerical  values  accordingly  are  extrapolations  albeit 
rather  safe  ones  because  T and  U are  approaching  constants  at  the  end 
of  the  calculation. 


The  ratio 
rB 


UBa° 


Uaf 


1/3 


= 5/f. 


1/3 


5.62 


(152) 


according  to  Equations  150  and  151.  This  implies  f = .70  at  the  end 
of  the  calculation  (since  r~5Ua)  and  therefore  the  final,  time  similar 
form  may  have  been  approached.  We  cannot  make  a stronger  statement  be- 
cause the  ratio  T/Ua  is  not,  at  present,  known  to  sufficient  accuracy. 


If  buoyancy  were  the  only  source  of  early  circulation  we  could 
write  a formula  for  X at  this  point.  However,  it  must  be  true  that 
blast  reaction  and  ground  interaction  interfere  with  the  mechanism  of 
circulation  generation  by  buoyant  forces. 


To  make  this  apparent,  and  to  find  a way  to  estimate  this  inter- 
ference, we  will  develop  a simplified  theory  of  buoyant  circulation. 


A simple  approximation  is  to  regard  the  vortex  as  spherical 
with  constant  radius,  a0 , until  it  rises  through  a distance  like  a0 
(i.e.,  until  it  turns  itself  through  about  tt  radians). 


The  impulse  equation  needs  to  represent  contributions  from 
fireball  materia]  as  well  as  the  exterior  flow  field, 

A [(p  + pa/2)a^  U]  = gaj  A . (153) 
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We  can  treat  everything  except  U and  t as  constant  in 
Equation  153,  then 


U = t /T,  t + U 

(p  + pa/2)  o 


(154) 


where  we  have  introduced  an  initial  value,  , to  represent  vortex 
velocity  imparted  by  some  other  mechanism,  such  as  blast  reaction  cr 
ground  interaction. 

The  rise  height  is  obtained  by  integrating  Equation  154, 


(Pa  + 2p) 


(155) 


The  torus  presumably  forms  and  cuts  off  generation  cf  circula- 
tion when  z~a0.  The  torus  time,  t,  is  obtained  by  equating  z to  aQ 
in  Equation  155  and  solving  for  t, 


(Pa  + 2^)Uo 


I 48aoA 
1 + (Pa  + 2p^  “ 1 


Limiting  behavior  of  Equation  156  is. 


lim  x = 


'(Pa  + 2p)a0 


(156) 


(157) 


lim  t = ao/U 


U -+  00 
o 


(158) 


The  formula  for  U at  torus  time  is 


U = t * ■ t + U 

(P  + Pa/2)  o 


(159) 


1 


and  the  corresponding  value  of  X is 


P‘/3  ga0A 

X = p,/3  a U = 2 

a o (p„  + 2p) 


T + Xr  + XG 


We  can  normalize  our  buoyant  contribution  to  X by  forcing  the 
limit  as  U ->■  0 to  match  the  results  of  Reference  6,  (Equation  150). 


Let 


P.1/3  ga  A 

XD  = N 

B (Pa  + 2p) 


(161) 


and  require 


lim  X = p‘<>  r /s 

■i-o  0 

O 


(162) 


Substitution  of  Equation  157  into  161  and  Equation  150  into 


162  vields 


\ = 4.5/5  = 0.9 


(165) 


after  we  take  the  limit  p -*■  0,  consistent  with  Equation  150. 


Our  formula  now  reads 


X = 


XB  + Xr  + XG 


(164) 


■.here 


X r,  = 0.9 


0l/3 

a 


;>:r  ganA 


(P_  + 2p) 

« 


(165) 
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and  T is  given  by  Equation  (156) 


The  effect  of  this  treatment  is  to  reduce  the  buoyant  contri- 
bution to  X to  some  extent  whenever  some  other  source  of  circulation 
acts  promptly,  the  tendency  will  be  to  partially  counterbalance  effects 
of  other  sources  on  the  long  term  rise  history. 

We  should  point  out  that  an  attempt  at  precise  treatment  of  the 
interaction  of  these  several  sources  of  circulation  is  formidable.  For- 
tunately, at  least  the  major  approximation  made  here,  that  both  the  blast 
reaction  and  ground  shock  occur  prior  to  the  time  significant  circulation 
has  been  buoyantly  generated,  is  a good  one  whenever  these  effects  are 
significant. 

10.0  ENDING  THE  CALCULATION 

One  has  the  intuition  that  the  time  similarity  assumption 
must  break  down  after  the  mass  discrepancy,  Aa3,  goes  negative  (i.e., 
when  the  fireball  becomes  more  dense  than  its  surroundings).  In  addi- 
tion, zero  buoyancy  is  a basic  condition  for  stabilization. 

One's  first  inclination  is  therefore  to  stop  calculating  when 
Aa3  equals  zero  and  conclude  that  the  correspondirg  altitude  is  approx- 
imately equal  to  stabilization  altitude. 

However,  if  one  looks  at  numerical  results  he  finds  that  the 
momentum  of  the  vortex  is  generally  adequate  to  carry  thr  vortex  well 
past  the  altitude  where  the  mass  discrepancy  goes  to  zero,  with  any 
rational  entrainment  assumption.  Further,  the  combination  of  entrain- 
ment and  expansion  reduce  the  vortex  density  with  adequate  efficiency 
that  it  cannot  fall  back  a very  large  fraction  of  the  overshoot  distance. 
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Rather  than  embark  on  an  effort  to  develop  a new  theory  for 
this  regime  we  choose,  for  the  first  attempt  at  least,  to  see  how  far 
the  present  theory  can  be  pushed.  Accordingly  we  will  continue  calcu- 
lating until  either  the  mass  goes  to  zero,  due  to  a negative  a,  or 
the  velocity  goes  to  zero  (X  = 0) , Stabilization  will  be  assumed  to 
occur  when  the  mass  has  been  reduced  to  half  its  original  value,  since 
we  are  after  mean  altitude,  or  when  X goes  to  zero  in  cases  where  a 
does  not  go  negative.  F;or  consistency  we  should  reduce  X as  mass  is 
detrained,  we  will  do  this  on  the  basis  of  cross  sectional  area  (a2). 
Perhaps  we  should  also  consider  the  effect  of  this  contribution  to 
dX/dz  on  a but  it  is  not  worth  the  effort  until  one  has  some  evidence 
that  the  model  resembles  reality  in  this  regime.  If  this  method  were 
correct  we  could  hope  to  find  not  only  the  debris  altitude  and  extent 
versus  time  during  the  rise,  but  also  the  final  vertical  extent  of  the 
stabilized  cloud.  The  lowest  point  would  be  where  a goes  negative  and 
the  highest  extent  where  the  mass  goes  to  zero.  One  would  be  very  sur- 
prised if  nature  were  so  cooperative. 
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